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Abstract 

We suggest a simple model, based on the type-I seesaw mechanism, for the lepton 
mass matrices. The model hinges on an Abelian symmetry which leads to mass 
matrices with some vanishing matrix elements. The model predicts one massless 
neutrino and = (M is the effective light-neutrino Majorana mass matrix). 
We show that these predictions perfectly agree with the present experimental data if 
the neutrino mass spectrum is inverted, i.e. if ma = 0, provided the Dirac phase 5 is 
very close to maximal (ib7r/2). In the case of a normal neutrino mass spectrum, i.e. 
when mi = 0, the agreement of our model with the data is less than optimal — the 
reactor mixing angle ^13 is too small in our model. Minimal leptogenesis is not an 
option in our model due to the vanishing elements in the Yukawa-coupling matrices. 
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Notation: In the model of this Letter we utilize the type-I seesaw mechanism. There 
are right-handed neutrinos z/r which are invariant under the SU{2) x f/(l) gauge group. 
One defines the left-handed conjugates z/^ of the through 

^'l = Cul (1) 

where C is the Dirac-Pauli charge-conjugation matrix. The neutrino mass terms are given 
by 

C,m-^ss = -URMDi^L-i^LMjji^R-^URMRCul+^ulC-^M-^UR (2) 

where, without loss of generality, the matrix Mr is symmetric. Under the conditions of 
the seesaw mechanism, the approximate Majorana mass matrix for the light neutrinos is 
then 

M = -MIM^^Md. (4) 
The symmetric matrix M is bi-diagonalized by a unitary matrix U via 

U'^MU = diag(/xi, /i2, /is), (5) 

where Hj = rrije^'^^^ (j = 1,2,3, no sum over j); the (real and non-negative) rrij are the 
neutrino masses and the /3j are their respective Majorana phases. The matrix U may be 
parametrized 

(C12C13 S12C13 sise"*"^ \ 

-S12C23 - Ci2S23Sl3e*'^ C12C23 - Si2S23Sl3e*'^ S23C13 , 
S12S23 - Ci2C23Si3e'^ -C12S23 - Si2C23Si3e*'^ C23C13 / 

(6) 

where the 7^ (a = e, /i,r) are unphysical phases, q = cos6'j {i = 12,23,13), Si = sin^j, 
and S is the physical (observable) 'Dirac' phase. The angles 6i belong to the first quadrant. 



Purpose: Recently [T], one of us has suggested a model that yields, in a certain limit 
and in the basis where the charged-lepton mass matrix Mg is diagonal, the predictions 

detM = 0, (7a) 
= 0. (7b) 

It has been claimed [1] that the predictions ([7]) "are compatible with the experimental 
data on neutrino masses and on lepton mixing" but this assertion has not been proved 
and qualified in ref. [1]. It is the aim of this Letter to 

1. Present a model, simpler than the one in ref. PT], that leads to the predictions 

2. Work out analytically the consequences of those predictions for the parameters of the 
matrix U, both in the cases of an inverted and of a normal neutrino mass spectrum. 
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3. Display graphically those predictions by means of scatter plots. 

4. Analyse the consequences of our model for leptogenesis. 

5. Show that, in the context of our model, the conditions ([7]) are stable under the 
(one-loop) renormalization-group evolution of M. 

To be sure, the conditions ([7]) have been studied in the literature before, see for instance 
refs. [21 [3]. However, the points 1, 4, and 5 above are original, and we believe that we 
have also performed the tasks 2 and 3 in a simpler, more transparent fashion than in the 
previous literature. 



Mass matrices: We shall assume that there are only two right-handed neutrinos vr [2]. 
This means that the matrix Md is 2x3 and, therefore, the 5x5 mass matrix in equation ([3]) 
has a 3 X 3 null submatrix, hence it has (at least) one zero eigenvalue, i.e. there is one 
massless neutrino. Furthermore, in the basis where Mr is diagonal, 

Mr = diag (x, y) , (8) 

the matrix has two texture zeros: 

By utilizing equation (jl]), one then has 

(a'^/x ac/x \ 

b'^/y hd/y . (10) 

ac/x hd/y c^/x + d'^/y j 

This matrix features equations (JTj). Notice that we are implicitly assuming that the 
charged-lepton mass matrix is diagonal too. The condition (17bl) means that 



/iIf/llf/21 + /i;f/l2f/22 + ^HUl^U2Z = 0. (11) 

The condition (ITall means that one of the neutrinos is massless. Thus, either fii = and 
then the neutrino mass spectrum is normal, with m2 = y^Am^ and = a/ Am^tm' 
or m3 = and then the neutrino mass spectrum is inverted, with mi = a/ Am^^^ and 
m2 = v^Aml~+Ar 



Inverted spectrum: In the inverted-spectrum case, equation ffTTl) implied 



'^lC?2 (^12^23 + Cusl^sj^ + 2C12S12C23S23S13 COS 6) 

^12 ('^12'^23 + '^12'^23'^13 ~ 2C12S12C23S23S13 COS 5) . (12) 



^The relative phases of the terms in the left-hand side of equation ^Tl\ depend on the Majorana phases 
and are therefore of little relevance. 
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Therefore, 



n _ ( 2 _ 2\ 2 2 2_|_/24_ 24\2 2 
U — ^l) '^12'^12'^23 ~r l^2'^12 ^l'^12) ■^23'^13 

- {mlsl2 + mlcl^) (2C12S12C23S23S13 COS 5) . (13) 

One may compute the approximate value of each term in the right-hand side of equa- 
tion ( fT3l) by using as a guide the Ansatz of tri-bimaximal mixing |1] S12 = 3"^^^, S23 = 
2-1/2 'jj^^g^ equation (fT3|) may be approximated by 



Ami, Am2^ , 
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^ ^i:^ - ^^a^ - ^ ^^""^'^ s,s COS 5. (14) 
g g 13 3 13 V ; 

Numerically, using 5^3 ~ 0.017 [5], the first term in the right-hand side of equation (IHj) 
is about 8.5 x 10~^ eV^, the second term is about —6.5 x 10~^ eV^, and the coefficient of 
cos 5 in the third term is about 1.4 x 10~^ eV^. One concludes from this estimate that 
equation (fT4l) imphes 

cos (5 ^ 0. (15) 

Remarkably, equation (JT3l) does not imply any constraint on the mixing angles — it works 
perfectly well for any of their phenomenologically allowed values. One concludes that, in 
the case of an inverted neutrino mass spectrum, our conditions ([7]) fit perfectly well the 
phenomenological values of the mixing angles and make the prediction f|T5|) for the Dirac 
phase. 

Normal spectrum: In the normal-spectrum case, equation ( fTTI) implies 

^"^soi^L {cucls + sL^L^L - 2C12S12C23S23S13 COS 6) = Aml^^sj^sl^. (16) 

In the left-hand side of this equation, the first term inside the parenthesis dominates over 
the other two because S23 ~ C23, cu > su, and S13 ^ 1. Therefore, 

^^sol ^ ^ ■^13 i^-^Y) 

Am^tm C23 C12S12 

Phenomenologically, the left-hand side of equation (fT7|) is about 0.03, while in the right- 
hand side S23/C23 ~ 1 and 1 /(cf2"5i2) ~ 4.5. Therefore, the model predicts 5^3 ~ 0.0067, 
which is somewhat below the best-fit value for 5^3 [5j. Moreover, while 1/(0^2-312) is 
rather stable over the 6^12 allowed range, sl^/c^^ lies at 3 cr in between 0.64 and 1.79, 
which means that sf^ must be significantly correlated with 533. One concludes that, if the 
neutrino mass spectrum is normal, then the condition fl7bp does not give a perfect fit to 
the data because it implies a much too low 613] moreover, there should be a correlation 
between the mixing angles 613 and 6'23. The phase 6 remains, in the case of a normal 
neutrino mass spectrum, free. 

Neutrinoless 2/3 decay: One is particularly interested in the quantity m^^ = |Mee|, 
which is the mass relevant for neutrinoless double-beta decay. One easily finds that 
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equation (fTTi) together with the unitarity of U imply 



ms = ^ iTLpp = nil 



mi = =^ misi3 = m^ 



U13U31 
U22 

U22 



~ mi, 




(18) 



(19) 



where in the second step we have taken into account the approximate values of the ma- 
trix elements of U given by the Ansatz of tri-bimaximal mixing. One sees that mpp ~ 



0.005 eV in the normal case [B]. Thus, in our model m^^ might be measurable if the 
neutrino mass spectrum is inverted, but is certainly too low to be measured (at least any 
time soon) if the mass spectrum is normal. 

Numerical analysis: We now undertake a scan of the parameter space allowed by the 
most recent phenomenological data on neutrino masses and oscillations [S]. As explained 
earlier, in the inverted-hierarchy case the fit of our model to the data is effortless and 
we only had to use for the observables the la intervals in ref. [5]. The fitting procedure 
consisted in varying the several observables (the angles 6*12, 6*13, and 6*23, the mass pa- 
rameters mi and m2, and the phase 5) over their allowed ranges, but making them obey 
equation 0131) . Remarkably, this fit to the experimental data enforces no constraints on 
the parameters of the matrix U , with the exception of the phase 5, which is found to be 
close to its maximal value ±vr/2, as seen in figure [H In that figure, we plot the quantity 
mpp = |Mee|, relevant for neutrinoless double-beta decay, versus the sine of 6. Each point 
in that plot corresponds to a particular combination of ^/-matrix parameters which obey 
the constraint of equation f lT^ . We thus confirm the two estimates shown previously: 
that the CP-violating phase 6 is close to its maximal value, and that m^^ ~ 0.05 eV. 

The fitting procedure is analogous for the normal-hierarchy case, except that we now 
want the parameters to obey equation ( IT6!) instead of equation ( [T3|) . It turns out that the 
fit is impossible if one uses the la intervals for the various parameters; we have therefore 
used the 3a intervals instead. The main problem is with the angle 613, which is driven by 
the condition (ITB|) to values much too small. This can be appreciated in figure [2], where 
we have plotted sin^ 6'i3 against sin^^23- The red cross indicates the phenomenological 
best-fit value; one sees that our model yields 61^ always well below that best-fit value. In 
figure [2] one can also observe the correlation between the ^13 and 623 mentioned earlier. 
We have found no other correlations between observables in this normal-hierarchy case. 
Unlike in the inverted-hierarchy case, the value of m^^ is found to be always extremely 
small — in between 0.002 eV and 0.0035 eV — and the CP-violating phase S can take any 
value. 

A model: We now explicitly construct a model which generates the neutrino mass 
matrices in equations ([8]) and ([9]) together with a diagonal charged-lepton mass matrix. 
Those mass matrices are solely characterized by some zero matrix elements, and it is 
known [7j that any set of texture zeros can always be justified through a model furnished 
with an Abelian flavour symmetry and, possibly, a large number of scalar gauge multiplets. 



'atm 



0.05 eV in the inverted-spectrum 
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Figure 1: Scatter plot of the Dirac phase 5 against m^^ = |Mee| in the case of an inverted 
hierarchy of neutrino masses. The fit is made to the la data of ref. [5]. Notice that the 
fit strongly pushes 6 to its maximal values. 

Indeed, there are always many such possible models; the model that we shall present in 
this paragraph is only one (possibly minimal in terms of its scalar content) out of many 
others which might perform the job. Our model is based on the Abelian symmetry Zg; let 
a = = exp {in/A). The model has three scalar gauge-5'[/(2) doublets $6 = (0^, 0°) 
{b = 1,2,3) and one complex gauge-invariant scalar S. Under Zg, $1 is invariant and 

Therefore, the scalar potential will contain the phase-sensitive terms ${$2 ^\^2, ^\^3 S, 
and S"^, together with their Hermitian conjugates, thereby avoiding the appearance of any 
Goldstone bosons when S and the 0° acquire vacuum expectation values (VEVs). The 
matrix Mr is generated by the Yukawa couplings of the right-handed neutrinos to S, 

^SYuk = -y uriCur^S - y iyR2CuR2S* + H.C., (21) 

yielding x = yis and y = y2S*, where s = (0 l^l 0) is the VEV of S. The left-handed lepton 
doublets Dia = {^La, oil)^ and the right-handed charged-lepton singlets aR transform 
under as 

DLe ctDlc, Dl^ ct^Dl^,, D lt cf'^Dlt, eR acR, fj,R a"^ fiR, tr -> tr. (22) 
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Figure 2: Scatter plot of the reactor angle against the atmospheric angle 623- The fit 
is made to the 3a phenomenological observables of ref. 0; the cross indicates the best-fit 
point found in that paper. Notice that our fit pushes 6*13 to smaller values. 

This leads to the following Yukawa Lagrangian: 

^Cd-Yuk = - (0?, -(Pt) {y^^RlDhe + yiVR2DLT) - (03, "03 ) {V^oVriDlt + 2/6'^H2^L/.) 

-yi {<Pi^ 4>T) ^rDlc - Vs (03 ' 03*) P'RDL^^ - 2/9 (02^, 0°*) t^rDlt + H.c. (23) 

Let Vb = (O|0|^|O) denote the VEV of 0°, then the charged-lepton masses are rrie = 
lyjVil, rrin = lysv^l, and rrir = 1 2/9^2]; the elements of the neutrino Dirac mass matrix 
Mo in equation are a = y-sVi, b = y^v^, c = y^v^, and d = y^Vi. Notice that, as 
mentioned earlier, the charged-lepton mass matrix is automatically diagonal with this 
symmetry. Furthermore, unlike what usually happens in flavour models based on non- 
Abelian symmetries (for instance, a recent proposal that we have made in ref. ^), there 
is here no specific requirement (beyond conservation of electromagnetism) on the model's 
vacuum and, therefore, the fit to the leptonic sector imposes no constraints on the scalar 
sector. 

Leptogenesis: We next consider leptogenesis. For reviews of leptogenesis see, for in- 
stance, [01 dn]. We firstly do this in the context of a general multi-Higgs-doublet model 
(MHDM). Let there be uh Higgs doublets (6 = 1,2, ... ,nH) and n^, right-handed 
neutrinos {i = 1,2, . . . , riy); the latter are supposed to be eigenstates of mass and, in 



7 



this paragraph, rrii denotes the mass of ujn. The Yukawa couphngs of the right-handed 
neutrinos are given by 



6=1 i=l a=e,fj,,T 



+ H.C., 



(24) 



where the Yj, are dimensionless riiy x 3 Yukawa-couphng matrices. In the approximation 
where the right-handed neutrinos are (very) massive but both the charged leptons and 
the charged scalars 0^ are massless, the CP- violating asymmetry between the decay rates 
of Ujii to a'^cj)'^ and a;~0^, summed over the flavour a, is given by 



1 



EE^- 

c=l j=l 



YY^ 



-|-Im 

where Xj = / ml and [S] 

f{x) 
9{x) 



y,yA (y;y;t) 



(25a) 
(25b) 



1 (1 + x) In 



X 



1 + x 



X 



We now consider the particular case of our model in the previous paragraph. There, 



2/3 

7/4 



Y, 









2/6 



(26) 
(27) 

(28) 



It is easy to see that the matrices YiY^ and ^3^3^ are diagonal, and therefore do not 

contribute to line fl25bl) . while (^lY^^ = and ^Yay/j = 0, and therefore line fl25al) 

vanishes too. We conclude that leptogenesis (at least in its simplest version, where the 
flavour of the lepton a is summed over) is not viable in our model. Indeed, this is a 
general result, i.e. it holds because of the texture zeros in our mass matrices and not only 
in our specific model of the previous paragraph. Since M/j is diagonal, we have mi 
and 7712 = \y\- We know that 



\x\ 



riH 



M, 



D 



(29) 



b=l 



In order for the in equation to follow from an Abelian symmetry which enforces 
its particular texture zeros, the conditions 



(n)i, 

(n)2e 

(n)i.(n)2. 



0, 
0, 




(30a) 
(30b) 
(30c) 
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[no sum over b in equation fl30cp ] must hold for all 6 = 1,2,..., uh- [If equation fl30cl) did 
not hold, that would mean that u^n and z//j2 transform in the same way under the Abelian 
symmetry, and then {Md)^^ = would imply {Md)^^ = 0, etc.] Now, equations (l30l) 
imply that 

^ Wi«(Vk = 0, (31) 
and therefore line (125bl) vanishes. Furthermore, because of equations (I30ap and (130b|) . 



(nr;)^, = (n)i^ (f;)2, (f.)i, (n*)2. , (32) 

which vanishes because of equation (I30cp . Therefore, line (I25ap is zero and = 0. 



Renormalization-group invariance: We now show that the conditions ([7]) are, in the 
context of the model that we have suggested in this Letter, invariant under the (one-loop) 
renormalization-group (RG) evolution of M. In that evolution, one must consider in 
the context of a MHDM, the effective operators 

0,c= E 4^)^-1 - (33) 

a,l3=e,ii,T 

where the k*^^^^ are matrices in flavour space. The effective operators Obc arise in the 
seesaw mechanism upon integrating out the heavy right-handed neutrinos. One has /t^^^-* = 
—Yf^ M^^Yc] in our model, from equations ([8]) and fl28l) . all the k^^^^ have, at high energy, 
vanishing (e, /i) and (/i, e) matrix elements. During the RG evolution to lower energies, 
the operators Obc mix among themselves. Moreover, the RG derivative dn^^'^^ / d(ln/i) 
contains [11] matrices n^'"^^ multiplied on the right by matrices 



"fee 



Y^t (34) 



or multiplied on the left by matrices F^-^. In equation the Yukawa-coupling matrices 
Yb are those which appear in the Yukawa couplings of the right-handed charged leptons, 

/:£Yuk = -V V (Yb) «R (0,-z/i/3 + 0r/3L) +H.C. (35) 

6=1 Q,f5=e,fi,T 

In our model the matrices — cf. equation (123!) — 

Yi = diag (1/7, 0, 0) , Y2 = diag (0, 0, y^, ) , % = diag (0, ys, 0) (36) 

are diagonal and, therefore, the Fbc are diagonal too. Hence, the multiplication of the k^'"^^ 
by the Fbc preserves the texture zeros that exist in the k^''^\ Since in our model all the k^'"^^ 
initially have null (e, /i) and (/x.e) matrix elements, that feature will be preserved all along 
the RG evolution. Therefore, the condition (l7bl) is RG-stable. As for the condition (|7all . it 
follows directly from the fact that there are only two right-handed neutrinos, and therefore 
it is stable independently of the specific form of the RG equations. 
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Conclusions: In this Letter we liave sliown tliat tlie conditions ([7]) on the effective hght- 
neutrino Majorana mass matrix M may be imposed by means of an Abehan symmetry 
in a three-Higgs-doublet model with an extra scalar singlet. We have furthermore 
demonstrated that those conditions agree perfectly well with the known phenomenology 
if the neutrino mass spectrum is inverted, in which case the Dirac phase 5 is predicted 
to be extremely close to ±7r/2. The conditions ([7]) do not work as well if the neutrino 
mass spectrum is normal, since in that case they necessitate a rather small ^13, but they 
cannot yet be excluded at more than the la level. Unfortunately, (flavour-independent) 
leptogenesis cannot work in the scenario in which the conditions ([7]) are enforced by means 
of an Abelian horizontal symmetry. 



Addendum: After completion of this work, our attention has been called to a recent 
paper [12j which has suggested, albeit from a different starting point, the conditions ([7j). 
That paper states that the condition (ITbl) "suffer [s] from little predictivity" in the case of 
an inverted neutrino mass spectrum; on the other hand, that case is, in our view, the one 
in which the condition fl7bl) agrees best with experiment. 
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